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ABSTRACT 


The Integrated Systems Program is capable of simulating both nominal 
and anomalous operation of the Apollo Cryogenics Storage System (CSS). Two 
versions-of the program exist; one for the Apollo 14 configuration and 
the other for "J" Type Mission configurations. The program consists of two 
mathematical models which are dynamically coupled. A model of the CSS com- 
ponents and lines determines the oxygen and hydrogen flowrate from each storage 
tank given the tank pressures and temperatures, and the Electrical Power Sub- 
system (EPS) and Environmental Control Subsystem (ECS) flow demands. Tempera- 
tures and pressures throughout the components and lines are also determined. 

A model of the CSS tankage determines the pressure and temperatures in the 
tanks given the flowrate from each tank and the thermal environment. The 
model accounts for tank stretch and includes simplified oxygen tank heater and 
stratification routines. Supporting subroutines are used to calculate thermo- 
physical properties from National Bureau of Standards (KBS) data and to calcu- 
late the CSS thermal environment during various mission phases. Anomalous 
conditions such as a fan and/or heater failure, degraded tank vacuum annulus 
performance, failed check values, and leaking lines can be simulated. The 
program is currently operational on the Univac 1108 computer. 
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1.0 INTRODUCTION 


The objectives of Task 705-2, "Apollo Service Module Integrated Cryo- 
genic Storage System Analysis," are to develop the analytical capability 
required to: 

a) _ Develop an integrated mathematical model of the Apollo Cryogenic 

Storage System (CSS) , 

b) Provide support to the redesigned CSS oxygen tank development 
program. 

c) Predict flight performance of the CSS to verify its adequacy to 
meet mission requirements. 

d) Carry-out postflight analysis of the CSS to confirm the adequacy 
of the mathematical model. 

To accomplish these objectives, the task was divided into six subtasks: 

Subtask 1 - Basic Cryogenic Tank Program 

Subtask 2 - Advanced Cryogenic Tank Performance Program 

Subtask 3 - Integrated System Program 

Subtask 4 - CSS Oxygen Tank Development Support 

Subtask 5 - Preflight Mission Analysis 

Subtask 6 - Plight Analysis 

This document describes the Integrated Systems Program developed under 
Sub task 3. 

The Integrated Systems Program is a mathematical representation of the 
Apollo Cryogenic Storage System (CSS) and is capable of simulating both 
nominal and anomalous operation of the system. The program is modular in 
structure and consists of a storage tank model, a components and lines model 
and auxiliary routines. Two versions of the program exist; one for the 
Apollo 14 spacecraft configuration and the other for the "J" Type Mission 
configuration. 
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The "J” conf iguratxon dxffers from the Apollo 14 configuration prxmarxly by 
the addition of a hydrogen storage tank and several other components. 

The tankage model xs capable of sxmulating the thermodynamxc performance 
of the CSS hydrogen and oxygen tanks. The model accounts for the effects of 
heat leakage, heater and/or fan operatxon, duty cycles, and control equxpment 
as well as all signxfxcant thermal, elastxc, thermoelastxc and thermodynamxc 
effects which xnfluence tank performance. Capabilxty to predxct the effects 
of thermal stratxf xcatxon on the oxygen tank performance and to define oxygen 
tank heater temperatures are included. 

The components and Ixnes model determxnes the oxygen and hydrogen flow- 
rate from each storage tank based on the pressures and temperatures of the 
indivxdual tanks and the flow demands of the Electrical Power Subsystem (EPS) 
and the Envxronmental Control Subsystem (ECS). A transxent analysis xs per- 
formed to determine the pressures and temperatures at each junction and the 
flowrates in each leg of the system. These values are used to determxne the 
instantaneous steady-state pressure drop across all components in each leg. 

Anomalous operatxon of the tanks, such as heater and/or fan faxlures, 
degraded tank vacuum annulus performance, and loss of tanks, may be simulated. 

In addition, line leaks and anomalous operation of any component may be 
simulated for the oxygen and/or hydrogen distribution systems. 

The auxxlxary routines descrxbe the thermal environment of the CSS and 
provide oxygen and hydrogen thermophysical data. The vehicle thermal model 
computes the outer wall temperatures of the storage tanks as well as space- 
craft bay wall and shelf temperatures for various mxssion phases Includxng 
translunar coast, transearth coast, attitude holds, and lunar orbit. The 
thermophysical data routine contains data for both oxygen and hydrogen ob- 
taxned from the National Bureau of Standards (NBS) over the expected system 
Operating range. 
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Input to the Integrated Systems Program includes the initial state and 

> « 

quantities of the cryogens in the tanks, EPS and ECS flowrate demands, com- 
ponent characteristics, inertial acceleratxons and the mission mode. 

The program, written in Fortran V, is operational on the Univac 1108 
computer^ 
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2.0 SYSTEM DESCRIFTION 


The Apollo CSS supplies oxygen and hydrogen on a demand basis to the EPS 
and the ECS. The cryogenic fluids are stored in their supercritical states 
in tanks located in the Service Module (SM) . Each tank is equipped with indi- 
vidual 'components for filling; venting; quantity, temperature, and pressure 
measurements; pressure control and relief; and the accompanying displays and 
controls. For all missions through Apollo 13, the storage system was comprised 
of two hydrogen tanks and oxygen tanks. The arrangement of these tanks in 
Sector IV of the SM is shown in Figure 2.1. For the Apollo 14 Mission the 
storage system consisted of an additional- oxygen tank located in Sector I of 
the SM making a total of three oxygen tanks and two hydrogen tanks. Subsequent 
missions will have three hydrogen tanks as well as three oxygen tanks. The 
additional hydrogen tank will also be located in Sector I of the SM. Figure 
2.2 presents the tankage arrangement in Sector I. 

During nominal operation, the storage system is maintained at a rela- 
tively constant pressure with fluid temperature increasing during the mission. 
Pressure is maintained above the critical pressure in the tanks by the 
application, either manually or automatically, of electrical heater power, ' 
Prior to the Apollo 14 Mission, electrical heat was supplied to the storage 
tanks in conjunction with destratification fans. However, subsequent to the 
Apollo 13 oxygen system failure, the system was redesigned to eliminate the 
fans in the oxygen tanks, although fans are still used in the hydrogen tanks. 

The nominal operating pressure of the storage tanks is 245 ± 15 psia 
and 900 ± 35 psia for the hydrogen and oxygen tanks, respectively. 

The fluid flow rate is controlled on a demand basis by the ECS and EPS 
regardless of the storage system's pressure or operating mode. Under nominal 
operating conditions there are 29.3 pounds of hydrogen and 330.1 pounds of 
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FIGURE 2.1 - SM SECTOR IV CSS TAMKAGE ARRANGE^^E^n 
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FIGURE 2.2 - SM SECTOR I CSS TANKAGE ARRANGEMENT 
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oxygen, loaded into each of the storage tanks. Of these amounts, there are 
28.2 pounds of usable 'hydrogen and 323.5 pounds of usable oxygen. Other 
pertinent storage system operational parameters are presented in Table 2.1. 

2.1 Cryogenic Storage Tanks 

A major difference between the oxygen and hydrogen tanks is in the 
material used in their fabrication. These materials are listed in Table 
2.1 and were obtained from Reference 1. In the broad sense, the physical 
configuration for both tanks is similar in nature and they can be described 
as a typical unit. 

Each tank consists of two concentric spherical shells with thermal 
insulation in the annular space between the shells. In the oxygen tanks, 
the insulation supports the pressure vessel. The annulus of each tank is 
equipped with a vac-ion pump to maintain the low pressure required for 
adequate vacuum insulation, although in flight the oxygen pumps are the only 
pumps which may be operated. In addition, the pressure vessel supports 
(hydrogen tank only), fluid lines, and the electrical conduit are located 
within the evacuated annulus. The exit fluid line is routed through the 
annulus before exiting from the tank -to provide vapor cooling. 

Prior to Apollo 13, each tank was equipped with a two element heater 
and two destratification fans. However, following the Apollo 13 oxygen 
system failure, the two element heater and destratification fans in the 
oxygen tanks were removed and replaced with a three element heater. In 
addition, the heater thermostat for protection against high heater tempera- 
tures was removed and replaced with a temperature sensor. The new sensor 
is located on the heater tube. Another temperature sensor is located above 
the capacitance probe and measures the bulk fluid temperature. Electrical 
lines leading into the redesigned oxygen tank have also been sheathed to 
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TABLE 2.1 


CRYOGENIC SYSTEM OPERATIONAL PARAMETERS 



Hydrogen 

Oxygen 

TANK MATERIAL 

5 Al-2.5 Sn ELI tl 

' Inconel 718 

TANK WEIGHT (PER TANK) 

- 


Empty (Approx) 

80.00 lb. 

90.82 lb. 

Usable Fluid 

28.15 lb. 

323.45 lb. 

Stored Fluid (100% 
indication) 

29.31 lb. 

330.1 lb. 

Residual 

4% 

2% 

Maximum Fill Quantity 

30.03 lb. 

337.9 lb. 

TANK VOLUME (PER TANK) 

6.80 ft^ 

4.75 

TANK PRESSURIZATION 



Heaters elements per tank 2 

3 

Flight Resistance 

Nominal Voltage 
Power 

78.4 ohms per 
element 
28 V DC 
10 watts per 
element* 

18.45 ohms per 
element 
28 V DC 

42.5 watts per 
element* 

Total Heater 
Heat Input Per tank 

68.2 BTU/Hr 

434.8 BTU/Hr 

Ground Resistance 

78.4 ohms per 
element 

18.45 ohms per 
element 

Voltage 

Power 

28 V DC 
10 watts per 
element* 

28 V DC 
42.5 watts per 
element* 

Total Heater 
Heat Input Per Tank 

68.2 BTU/Hr 

434.8 BTU/Hr 
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TABLE 2.1 (Continued) 



Hydrogen 

Oxygen 

Pressure Switch 



Open Pressure Max. 

260 psia 

935 psia 

Close Pressure Min. 

223 psia 

865 psia 

Deadband Min. 

10 psia 

30 psia 

Destratification Motors 
(2 Motors Per Tank) 



Voltage 

115/200 V 

None 


400 cps 


Power - Average 

3.5 watts per 
motor* 

None 

Total Average 
Motor 

Heat Input Per Tank 

23.8 BTD/Hr 

None 

SYSTEM PRESSURES 

— 


Normal Operating 

245 + 15 psia 

900 + 35 psia 

Spec Min. Dead Band 

10 psi 

30 psi 

of Pressure Switches 



Relief Valve Note: Relief Valves are referenced to environmental 

pressure, therefore pressure at sea level 

1 (psig) will be the same 

value in vacuum (psia) 

Crack Min. 

273 psig 

938 psig 

Full Flow Max. 

285 psig 

1010 psig 

Reseat Min, 

268 psig 

865 psig 

Outer Tank Shell 
Burst Disc 

+10 

■ 

Nominal Burst Pressure 

-20 psid 

75 + 7.5 psid 

SYSTEM TEMPERATURES 



Stored Fluid 

-425 to aO^F 

-300 to 80°F 

*Conversion Factor: 1 watt = 3.41 BTU/Hr 

» 
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reduce the fire hazard. 


2.2 Cryogenic System Components and Lines 

Figures 2.3 and 2.4 are schematics of the component and line configura- 
tions for the Apollo 14 oxygen and hydrogen systems, respectively, showing 
the relative locations of the components which make up the system. Figures 
2.5 and 2.6 present similar schematics for the *'J" Type Mission system. In 
general, the system is composed of lines, valve modules, flow filters and 
flow restrictors. For missions prior to Apollo 14, the system valve modules 
for both the hydrogen and oxygen systems each contained two relief valves, 
two pressure transducers, two pressure switches and one check valve. For 
Apollo 14, a half valve module was added for oxygen tank 3 and contains one 
relief valve, one pressure switch and a check valve. A solenoid operated 
- shutoff valve was also added between oxygen tanks 2 and 3 . 

The relief valves in these .modules are a gradual opening differential 
pressure poppet tjrpe which assure the lowest possible leakage with virtually 
zero delta pressure between crack and full flow. These valves are designed 
to be unaffected by back pressure in the downstream plumbing and are tempera- 
ture compensated over the full range of fluid temperatures which may occur. 

The pressure transducers are absolute pressure devices which utilize 
bonded strain gauges to produce millivolt signals to the signal conditioner. 
The signal conditioner is integral to the transducer with a 0 to 5 VDC out- 
put linearly proportioned to tank pressure. The transducers can operate at 
environmental temperatures as low as -80°F and are capable of supplying a 
continuous readout of pressure over the specified range with a maximum error 
of ± 2-1/2 percent of full scale. 

The pressure switches are single-throw, single pole absolute pressure 
devices which may be used to control the motor switches which activate the 
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destratification motors and/or tank heaters. These switches have been used 
■* « 

successfully in an environment of -400°F. 

The check valve is a spring loaded poppet type and is designed to open 
at a differential pressure of approximately 1 psi. The valve seat has a 
large area to prevent chattering during flow in the normal direction and to 
help obtain a positive seal if pressurized in the reverse direction. 

The fuel cell valve module for the hydrogen system consists of two check 
valves and three solenoid operated valves in one valve body. The check 
valves are a dual seat type which allow the main passage to be full open and 
the auxiliary passage barely cracked at low flows. At high flows, both the 
main and auxiliary passages are completely open. The solenoid valves are poppet 
latch type and are actuated by a magnetic armature suspended on a Bellville 
spring. Opening and closing coils on both sides of the armature actuate the 
valve. The valves open against pressure and pressure helps to seal the valve 
against leakage when closed to flow in the normal flow direction. 

For the oxygen system, the fuel cell valve module contains check valves 
similar to those in the hydrogen fuel cell valve module. However, subsequent 
to Apollo 13, the solenoid valves have been taken out of the module and 
mounted externally and upstream of the flow transducers to the fuel cells. 

The oxygen and hydrogen inline filters consist of several chemically 
etched discs stacked in a cartridge. They are rated at 5 microns nominally 
and 12 microns absolute with a capacity to contain .25 grams contaminant. 

Several components have been added to the CSS for missions subsequent 
to the Apollo 14 ("j" type missions). A half valve module containing a 
relief valve, a pressure switch, and a check valve has been added for the 
third hydrogen tank. There are three oxygen restrictors in the ECS for "j** 
type missions versus two for other Apollo missions. The third restrictor 
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line IS a branch of the oxygen tank 2 ECS supply line and reconnects down- 
stream of the junction of the other two restrictor lines. This restrictor 
line includes a filter and check valve. 


16 



3.0 ANALYTICAL FOKMULATION 

*• 

The Integrated Systems Program was developed as a tool to axd in the 
preflight planning and the postflight analysis of the Apollo Cryogenic 
Storage System. Modular in structure, the program is comprised of a tankage 
model, a components and line model and several auxiliary subroutines. All 
significant effects which influence system performance, such as thermal 
environment and mission mode, are Included. 

Given the component characteristics, the initial state of the cryogens 
in the tanks, the EPS and ECS flow requirements and the mission mode, the 
program is capable of simulating the nominal performance of the CSS. Figure 
3.1 presents a flow diagram of the program. Output includes the temperature, 
pressure and flowrates at various points in the system as well as the remain- 
ing quantities and state of the cryogens in the tanks. 

In addition to nominal performance, the program is capable of simulating 
the operation of the CSS under anomalous conditions. The effects of degraded 
performance of the tank vacuum annulus caused by the presence of air or the 
cryogen in the jacket may be determined. Abnormal environmental temperatures 
or tank fan and heater failures may be simulated. Anomalous operation of 
the check valves, filters, relief valves or flow restrictors may be investi- 
gated. Oxygen or hydrogen line leaks may be modeled nearly anywhere in the 
system. 

Two versions of the program were developed. One version is capable of 
simulating the CSS up to and including the Apollo 14 configuration. The other 
is capable of simulating the "J" T 3 q)e Mission configuration. The primary 
difference in these programs is in the components and lines model. Only the 
"J" Type configuration will be discussed in this document. 
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FIGURE 3.1 


INTEGRATED SYSTEMS PROGRAM 
FLOW DIAGRAM 
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The followxng sections discuss the analytical formulation of the ma 3 or 
models and supporting subroutines which comprise the Integrated Systems Program. 

"4 

3.1 Storage Tank Model 

The storage tank model was developed to simulate the thermodynamic 
performance of the Apollo CSS oxygen and hydrogen tanks . Given the initial 
conditions and flow demands, the model simultaneously updates the thermodynamic 
state in each storage tank based on control equipment operation, tank character- 
istics and heat leak. A flow diagram of the model is presented in Figure 3.2. 

The model may be divided into two sections. The first section deter- 
mines the maximum time increment over which all tanks may be simulated . The 
second section uses this time increment to update the state of the fluid in 
each tank corresponding with the conditions which prevail at the end of the 
time increment. 

The time increment depends on the rate of change of pressure and the 
operation of the pressure switches in each tank. The time for the tank 
pressure to reach a specified excursion (currently 15 psi) from its initial 
state IS compared to the time required to reach the pressure control limit. 

The shorter of these two time intervals is defined as the maximum dynamic 
simulation time increment for a given tank. Dynamic simulation time incre- 
ments for all tanks are compared and the smallest one is chosen as the total 
tank system dynamic simulation time increment. This value is then compared 
with a user specified increment and the lesser is selected as the time 
increment for use in updating the state of all cryogenic tanks . 

Upon the determination of the time increment, the thermodynamic state 
of the fluid in each tank is updated. The average pressure for the time 
interval is determined from the initial tank pressure, the initial pressure 
rise (or decay) rate, and the length of time interval. The average density 
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FIGURE 3.2 


STORAGE TANK MODEL 
FLOW DIAGRAM 
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FIGURE 3.2 


STORAGE TANK ilODEL 
FLOW DIAGRAJl 
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for the time interval is determined from the tank volume, the initial weight 
of cryogen in the tank, and the flowrate. Using these averaged thermodynamic 
properties, the average pressure change rate is recomputed for each tank. 

The average pressure change rate and flowrate are then used to update the 
pressure and density at the end of the time interval. Finally, fluid 
temperatures are obtained using the thermophysical routines. 

In addition to updating the state of each storage tank, the additional 
flowrates required by the fuel cells to provide power for all electrical 
equipment operating during the time interval are determined and added to 
the EPS flow demand. The time required for a complete pressure cycle and 
heater cycle is also computed for each tank. 

One of the basic functions of the model is the determination of the 
rate of change of tank pressure with respect to time. The method of deter- 
mining this derivative differs for the two cryogens since the effects of 
thermal stratification are considered significant in the oxygen tanks , but 
not in the hydrogen tanks. 

3.1.1 Oxygen Tank 

The thermal stratification, resulting from non-uniform heating of the 
cryogen, significantly affects the operation of the oxygen tanks. The non- 
uniform heating leads to higher pressurization rates than that which would 
result from uniform heating. Furthermore, at high densities, thermally 
stratified oxygen can experience rapid pressure decay upon mixing due to the 
thermodynamic properties of the oxygen. 

In order to determine the pressure derivative, the oxygen tank is 
represented by three nodes or volumes. One node is associated with 

the fluid within the tank heater column. A second node (V^) represents 
the fluid immediately surrounding the heater column. The final node (V^) 
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represents the remalnxng fluid in the tank. The nodal representation of 
the tank is shown in Figure 3.3. 

The oxygen which flows from the tank is withdrawn from Mass inter- 
change IS permitted between and but not between ^ ^ result 

of this constraint is that V. continuously increases during a mission due to 
fluid withdrawal from until a maneuver which would mix the fluid occurs. 

At that time the densities in the three volumes are assumed equal and the 
relative position of the boundary between volumes two and three are re-adjust' 
ed for this condition. The effects of tank stretch and the compression of 
the fluid in the fill and vent lines (see Appendix II) on the pressure 
derivative are accounted for in V^. 

The cylindrical heater is described by a heat source of a specified 
temperature with a specified heat flux into the two adjacent oxygen nodes. 


The heaters are mounted so that their axes are parallel to the CSM x-axes 
and perpendicular to the direction of the centripetal acceleration during 
passive thermal control rolls . Because of this , the heaters may be modeled 
as a horizontal cyclinder in a gravitational field. Based on this assump- 
tion, the net heating rate, q^, which is the sum of the convective and 


radiative heat transfer rates, may be derived as; 


q, = 525 irL. 
h h 


3,3 2 


k p gy c 




1/4 


Where g is the acceleration level and a is Stephan's radiation constant. 
This expression is based on a formulation presented in Reference 12. The 
heater temperature is then obtained by integrating the energy balance 
relationship : 
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FIGURE 3.3 

OXYGEN TANK NODAL MODEL 



I 


Heater 

Colmnn 



dt 

t 

This model of the oxygen tank heaters yielded excellent correlation with 
Apollo flight data. The heat flux to the surrounding nodes is calculated 
from the electrical energy input to the heater, the heater configuration, 
and the time history of the heater temperature Heat input is apportioned 
between the nodes adjacent to the heater by a simple conduction analogy 
which satisfies a total energy balance for the heater column. Heat conduc- 
tion IS permitted betxvfeen and , but not directly between and 
The thermal energy exchange between and must pass through the heater 
tube All external heat leak is applied to V^. 

Mass and energy balance relationships were developed for the fluid in 
each of the tank nodes. The energy balance equations are 

^ + ^1 [H (mj^)h2 + H (-m^)h^] (1) 

-m^ (h^ - 0^)} 



( 2 ) 


(3) 
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The mass balance equations are* 

*1 = “2 

^withdrawal 

The Nomenclature used in these equations is defined in Appendix I. 

The mass balance relations may be substituted into Equations (1) , 

(2) , and (3) to yield a set of three non-linear equations in the three 
dp . ^'^2 

variables lii^ and — ^ which may be solved iteratively. The resulting 
pressure derivative is then applied as described above. 

The pressure collapse potential is defined as the difference between 
the instantaneous value of pressure and the value that would result from 
sudden complete mixing of the fluid. Using the individual values of enthalpy 
for the three volumes, a mean enthalpy is determined for each tank. The 
mean density of the fluid in each tank is then determined from the known 
weight of the remaining oxygen and tank volume. The tank pressure for this 
mixed condition is obtained from the thermophysical routines. Finally, the 
pressure collapse potential is computed using this pressure and the pressure 
under stratified conditions. 
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3.1.2 Hydro Ren Tank 


During the Apollo 7 Mission, flight tests were performed which demon- 
strated that thermal stratification did not significantly affect the per- 
formance of the hydrogen tanks. It was therefore assumed that thermodynamic 
equilibrium existed in the hydrogen tanks at all times. The equation which 
determines the rate of change of pressure with respect to time for the 
hydrogen tanks is: 


iE. = 

dt 


64 . 

V 


m 


1 - 


3ap 


_3p_ 


V(1 + 3 gp I") 



6(|)p 9p 
l+3ap 8T 


"vl 3r 
VNp 2BY 


(1 - V) + 3a If 


(4) 


Equation (4) is derived in Appendix II and accounts for variations in tank 
volume due to fluid pressure and temperature excursions as well as the com- 
pression of fluid in the fill and vent lines . All terms of Equation (4) 
may be determined from the CSS system and fluid thermodynamic properties. 

3.2 Components and Lines Model 

The components and lines model simulates the operation of the Apollo 
CSS from the tanks to the user subsystems. Given the EPS and ECS flow 
demands, and the tank conditions, the pressures, temperatures and flowrates 
are determined at various points in the system. A flow diagram of the 
model is shown in Figure 3.4. 

The model is divided into two parts which are called, for convenience, 
the matrix and nodal models. The matrix model is a mathematical represen- 
tation of the oxygen and hydrogen systems based on the linearized electrical 
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FIGURE 3.4 


J TYPE COMPONENT AND LINES 
MODEL FLOW DIAGRAM 
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analogs shown in Figures 3.5 and 3.6, respectively. A transient, com- 
pressible flow analysis of each system is performed to determine the flo;? 
rate in each leg and the pressure at each matrix node. The nodal model 
then uses this information to perform a steady-state, incompressible flow 
analysis of each leg to determine the pressure at selected locations along 
that leg and to update the nonlinear characteristics of the individual com- 
ponents and line segment resistances. The nodal model also performs a 
transient thermal analysis of the lines and components. Detailed electrical 
analogs of the oxygen and hydrogen systems are shown in Figures 3.7 and 
3,8, respectively. Constants used in the system are presented in Tables 

3.1 and 3.2. 

3.2.1 Matrix Model 

The conservation of mass principle is used to determine the matrix 
node pressures. For any node i, this may be written (see Appendix III): 



dE. 


i dt 



(5) 


where : 


C 


i 


V 

1 

2RT^^ 


and; 



'iv®T 


RT. 

1 


dT. 

1 

dt 


( 6 ) 


(7) 


The matrix model solves this set of simultaneous equations for each system 
by an implicit backward differencing technique. The flow rate in each leg 
is then evaluated using the conductances in the characteristic matrix. 
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FIGURE 3.6 

ELECTRICAL ANALOG OF THE HYDROGEN SYSTEM FOR THE MATRIX MODEL 
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FIGURE 3.7 

THE OXYGEN SYSTEM FOR THE NODAL MODEL 
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TABLE 3.1 

J TYPE CONFIGURATION 
OXYGEN SYSTEM CONSTANTS* 


LEG 

NO. 

NODE 

NO. 

COMPONENT 

^EFF 

(in) 

D 

(in) 

2 

(Ibf/in ) 

Pr 

(Ibf /in^) 

SINK 

INDEX 

B 

mm 


m 

0.21 



151 

■ 


Filter 


0.21 







39.1 

0.21 



151 

■ 








2 

5 


36.0 

0.21 



151 


6 

Filter 

21.0 

0.21 





7 


37.1 

0.21 



151 


8 







3 

9 


18.3 

0.21 



163 


10 

Filter 

21.0 

0.21 





11 


42.3 

0.21 



163 


12 







4 

13 

Check Valve 

2.2x10® 

- 

3.0 

1.5 



14 







5 

15 



0.21 



163 


16 



0,21 



163 


17 



0.21 



159 


18 



0.21 



159 


19 


47.3 

0,21 



159 


20 







6 

21 

Check Valve 

1.4x10® 


3.0 

1.5 

151 


22 







D 

23 


40.0 

0.21 



151 

H 

24 


16.6 

0.21 



151 


25 


86.7 

0.21 



152 

■ 

26 


24.3 

0.21 



152 

■ 

27 


34.0 

0.21 



158 

■ 

28 








*Data obtained or derived from system drawings and References 2 and 3 . 
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TABLE 3.1 (Continued) 


LEG 

NO, 

node 

NO. 

COMPONENT 

^EFF 

(xn) 

D 

(in) 

Pc 

(Ibf/i-n^) 

2 

(Ibf/in ) 

SINK 

INDEX 

8 

29 

Shutoff Valve 

34.8 

0.21 



158 


30 



- 




9 

31 


21.2 




158 


32 


48.1 




153 


33 


17.3 

0.21 



160 


34 


58.8 

0.21 



160 


35 


13.8 

0.21 



160 


36 


79.3 

0.18 



161 


37 


63,4 

0,18 



161 


38 


28.4 

0.18 



161 


39 


24.2 

0.21 



162 

10 

41 


16.4 

0.21 



162 


42 

Filter 

21.0 

0.21 





43 


20.3 

0.21 



162 


44 

Restrictor 







45 


38.3 

0.21 



162 


46 







11 

47 

Check Valve 

3.2x10^ 

0.21 

3.0 

1.5 

162 


48 







12 

49 


28.5 

0.21 



162 


50 

Filter 

21.0 

0.21 





51 

Restrictor 







52 







13 

53 

Check Valve 

3.2x10^ 

0.21 

3,0 

1.5 

162 


54 
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TABLE 3.1 (Contxnued) 


LEG 

NO. 

NODE 

NO. 

COMPONENT 

^EFF 

(xn) 

D 

(xn) 

(Ibf/xn^) 

Pr 

(Ibf/xn^) 

SINK 

INDEX 

14 

55 


BHH|| 

mm 



162 


56 






162 


57 


77.9 

-0.21 



162 


58 

Shutoff Valve 






15 

59 


11.3 

0.21 



162 


60 







16 

61 


18.8 

0.21 



151 


62 







17 

63 


25.1 

0.21 



151 


64 


87.6 

0,21 



152 


65 


20.3 

0.21 



152 


66 


77.6 

0,21 



153 


67 


17.3 

0.21 



160 


68 


60.1 

0.21 



160 

- - 

69 


13.8 

0.21 



160 


70 


82.9 

0.18 



161 


71 


61.0 

0.18 



161 


72 


35,3 

0.18 



161 


73 


15.2 

0.21 



162 


74 

Filter 

21.0 

0.21 





75 

Restrictor 


- 





76 


6.9 

0.21 



162 


77 







18 

19 

Check Valve 

3.2x10^ 

0.21 



162 


WM 







19 

80 


15.6 

0.21 



162 


81 
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TABLE 3.1 (Continued) 


LEG 

NO. 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 


NODE 

NO. 

COMPONENT 

^EFF 

(in) 

D 

(xn) 

Pc 

(Ibf/in^) 

Pr 

(Ihf/in^) 

SINK 

INDEX 

82 

Check Valve 

3.2x10^ 

0.21 

3.0 

1.5 

162 

83 







84 

Shutoff Valve 

27.7 

0.21 



162 

85 







86 

Shutoff Valve 

50.1 

0.21 



162 

87 







88 


18.3 

0.21 



151 

89 

Filter 

21.0 

0.21 




90 


45.0 

0.21 



151 

91 







92 

Check Valve 

3.4x10^ 


3.0 

1.5 


93 



“ 

— 



94 

Filter 

21.6 

0.21 



152 

95 


24.9 

0.21 




96 







97 

Check Valve 

3,8x10^ 


3.0 

1.5 


98 







99 

Relief Valve 

48.1 

0.21 

999 

977 

151 

100 







101 

Relief Valve 

56.6 

0.21 

999 

977 

151 

102 







103 


43.6 

0.34 



151 

104 

Relief Vent 






105 

Relief Valve 

31.5 

0.21 

999 

977 


106 







107 


63.5 

0.21 



163 

108 

Relief Vent 






109 

Relief Valve 

49.8 

0.21 

1070 

1020 


110 


















































































TABLE 3,2 


J TYPE CONFIGURATION 
HYDROGEN SYSTEM CONSTANTS* 


LEG 

NO. 

■Hi 

B 

COMPONENT 


D 

(in) 

Pc 

(lb£/in^) 

Pr 

(Ibf /in^) 

SINK 

INDEX 

B 

■PM 


61.9 

0.21 



155 

B 

— 







2 

203 


74.5 

0.21 



154 


204 







3 

mfli 


59.5 

0.21 



163 


■■■1 







B 

207 

Check Valve 

4.0x10^ 


3.0 

1.5 


B 

208 







5 

209 


52.6 

lEI 



163 


210 


63.9 

wSm 



163 


211 


36.3 




159 


212 


51.8 

0.21 



159 


213 


70.1 

0.21 



159 


214 


58.2 

0.21 



152 


215 


44. 6 

0.21 



152 


216 


45.9 

0.21 



156 


217 


66.3 

0.21 



156 


218 







6 

219 

Check Valve 

4.0x10^ 

0,21 

3.0 

m 

156 


220 





■B 


B 

221 


8.3 

0.21 



156 

B 

222 

( 

22.9 

0.21 



155 


223 

Filter 

21.0 






224 


28.5 

0.21 



155 

B 

225. 







8 

226 

Check Valve 

126. xlO^ 


3.0 

1.5 



227 








*Data obtained or derived from system drawxngs and References 2 and 3. 
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TABLE 3 . 2 (Continued) 


LEG 

NO. 

NODE 

HO. 

COMPONENT 

^EFF 

(in) 

D 

(in) 

Pc 

(Itf /in^) 

Pr 

(Ibf/in^) 

SINK 

INDEX 

9 

228 



.0.21 



155 


229 

Filter 







230 



0.21 



155 


231 







10 

232 

Check Valve 



3.0 




233 





■H 


11 

234 

Relief Valve 


0.21 

280 

274 

155 


235 







12 

236 

Relief Valve 


0.21 

280 

274 

155 


237 







13 

238 



HSi 



— 


239 



Is&l 





240 






157 


241 






152 


242 



0.34 


- 

158 


243 



0,34 



158 


244 







14 

245 

Relief Valve 



280 

274 



246 







15 

247 



0.21 



BQI 


248 



0.21 





249 






,■11^ 

16 




0.34 



159 



Relief Vent 







37 


















































The program can. simulate up to two leaks in each system at any loca- 
tion except downstream of relief valves. Adding a leak is equivalent to 
adding a matrix node. These are included in the characteristic matrix as 
floating columns and rows. The flow at leaks is assumed to be choked and 
is evaluated in the same manner as overboard vents which are discussed 
below. 

3.2.2 Nodal Model 

The nodal model performs a steady state analysis of each leg using the 
matrix model node pressures and temperatures as boundary conditions. The 
individual components and line segments are analyzed as follows; 

3.2. 2.1 Lines and Fittings 

The pressure drop in a line segment from inlet i to outlet j is cal- 
culated by (see Appendix III) : 



128 LyR 
4 

™ 


(T^ + T ) m 


where: 


T. 

3 



UA 

me 

P 


+ T^exp 



( 8 ) 


(9) 


Since fluid velocities are small, it is assumed in the derivation of Equa- 
tion (8) that the flow is laminar and incompressible. The conservation of 
energy principle is used to determine the line temperatures. It is assumed 
that conduction through the stand off brackets and along the lines is neg- 
ligible. Then, for any node, k, 




k dt 




) - UA AT^ 


( 10 ) 
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Fittings have an equivalent L/D of approximately 3. Bends are included by 
an equivalent L/D which is a function of the bend angle (Figure 3.9). 
Clogged lines are simulated by incorporating an increased equivalent length 
into Equation (8) • 

3. 2. 2. 2- Filters , Check Valves , Relief Valves » and Shutoff Valves 

The pressure drops across filters and all open valves are evaluated 
from Equations (8) and (9). Based on typical values, an equivalent L/D 
of 100 is assumed for filters, open check valves, and open relief valves. 
The equivalent L/D for open shutoff valves is assumed to be 10. The pres- 
sure drop across all closed valves is calculated from: 


P “ P 


m 


,-7 


( 11 ) 
-10 


where the conductance G, is assumed to be 10 for check valves and 10 
for relief valves. Check valve and relief valve hysteresis is modeled as 
sho\<m in Figures 3.10 and 3.11, respectively- Clogged filters and valves 
are simulated by incorporating an equivalent L/D that is larger than the 
nominal value. Failed check valves and relief valves are simulated by in 
corporating a cracking and/or reseat pressure that will hold the valve 
open or closed as desired. 

The temperature of filters and individual valves are determined from 


, . dT 


-UA AT. 


LM 


( 12 ) 


while the temperatures of valves contained in a valve module housing are 
determined from 


(me ), 
p k dt 




(13) 
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FIGURE 3.9 

EQUIVALENT L/D FOR BEND 
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3. 2. 2. 3 Restrictors 


The pressure drop across the oxygen system flow restrictors Is evalu- 
ated using the data in Figure 3.12. It is assumed that the process Is 
Isothermal. 

3.2.2,4 “ Overboard Vents 

For compressible flow such as overboard venting the flow is choked. 
The pressure is evaluated by; 


P 


1 


v/*L 

K 


(14) 


1 

The nodal model also checks the solution for convergence. The updated 
resistance of each component and line segment is compared to the resistance 
used to calculate the flow rates. If the difference in resistor values is 
within the convergence criteria, the program proceeds to the next step time; 
if not, the flow rates are re-evaluated based on the updated resistances. 
Iteration of the characteristic matrix proceeds in this manner until the 
solution converges (normally within two iterations) . 


3.3 Auxiliary Routines 

In order to adequately simulate the operation of the Cryogenic Storage 
System, it was necessary to develop several auxiliary routines. Two of the 
most important routines will be discussed in this section. They are the 
vehicle thermal model and the thermophysical properties routine. 

3.3.1 Vehicle Thermal Model 

The vehicle thermal model was developed to compute the Cryogenic Stor- 
age System environmental temperatures and to determine the heat leak of the 
oxygen and hydrogen tanks under normal conditions or under conditions of 
loss of vacuum in the insulation annulus. The heat leak into the storage 
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tanks IS primarily a function of the flowrate through the vapor cooled 
shroud (VCS) and the temperature of the bay. The bay temperature depends 
to a large extent upon whether the vehicle is in an inertial hold mode or 
a passive thermal control mode. The effect of vacuum loss in the jacket is 
to increase the heat leak into the tank. 

3. 3.1.1 Oxygen Tank Heat Leak 

The heat leak into the oxygen tanks during nominal operating conditions 
is computed in the thermal model using empirical data correlations developed 
by the Beech Aircraft Company (Reference 4) . The correlations relate the 
heat transfer into the oxygen tank with the temperature of the oxygen in 
storage, the oxygen expulsion rate and the tank outer wall temperature. These 
correlations are given by the following equations : 

% - ALI _ A „ ( 15 ) 

- h'”') 

q = AL3 (T^ - T^) + .9 (17) 

The value of the constants which apply to Equations (15), (16), and (17) 
are as follows: 

Redesigned Tank: 


ALI = 

.9205 X lO”^ BTU/HR°R^’^ 

AL2 = 

.7278 X 10“^ BTU/HR°R^'^ 

AL3 = 

.03775 BTU/HR°R 


The amount and temperature of the oxygen in storage is computed in the 
tankage model of the Integrated Systems Program and is input to the thermal 
subroutine. Figure 3.13A presents typical oxygen tank heat leaks for various 
flowrates and quantities 
CHANGED JUNE 1972 44 



The outer wall tank temperature is computed within the thermal sub- 
routine through a series of data curve fit correlations. The data for 
these correlations were obtained using the Thermal Mathematical Model (TMM) 
of the Apollo Service Module. This program is composed of over 800 lumped 
parameters that constitute a thermal network specifically formated for 
solution with the SINDA Multi-Option computer program Using TMM, a set 
of parametric computer runs was performed varying the mission mode (inertial 
hold, passive thermal control or lunar orbit), the quantity of the cryogexi 
in the tanks, and the flowrates from the tanks A uniform internal tank 
temperature was assumed in determining the temperature of each oxygen and 
hydrogen tank surface. 

For use in the thermal subroutine, the temperatures computed for the 
surface of the oxygen tanks by the TMM were averaged for each tank. The 
average temperature were curve fit as functions of time. The form of the 


resulting curve fit equations is : 

2 3 

T = a + a-x + a„x + a_x 
o 1 2 3 


where x is the time argument. The coefficients for the polynomial curve 
fit equations are shown in Table 3.3. Using these equations, the tank 
outer wall temperatures are first estimated as functions of mission mode* 
time in mode and quantity. These values are then adjusted for the sun 
look angle. The temperature of the fluid in the VCS is then computed from 
the temperature of the fluid in storage and the outer wall temperature 
using an itei’ation scheme. Once the temperature of the VCS fluid has been 
determined, the heat flux is computed using Equations (15), (16), and (17) 
3 3.1 2 Hydrogen Tank Heat Leak 

The hydrogen tank heat leak for nominal conditions is computed 
in the thermal model using empirical data correlations developed by 
the Beech Aircraft Company (Reference 13) . These correlations are 
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TABLE 3,3 

CURVE FIT COEFFICIENrS REPRESENTATIVE 
OF THE CSS THERMAL ENVIROHMEKT 


Mission 

Mode 

Node 

a 

O 

CT) 

^1 

("F/llr) 

^2 

(°F/Hr^) 

^3 

(‘F/Ht^> 

T 

inf 

“F 

Time Araument 

t+3 

(Hrs) 

mm 

PTC 

Tanks 1 A 2 

A6 152 

35 715 

105 97 

-38 106 



\ 

PTC 

Tank 3 

25.164 

73 761 

202 92 

-112.08 



X 

Hot Soak 

Tanka 1 i 3 

60.0 

4 44 

,0 

0 0 


X 


Hot Soak 

Tank 2 

58 523 

31 715 

-3 57 

0 1277 


X 


Cold Soak 

Tanks 1 S 3 

-15 824 

368 40 

-485 5 

232 26 



X 

Cold Soak 

Tank. 2 

-61 612 

549 39 

-650 84 

283,34 



X 

PTC 

150 

66 0 

0 0 

0 0 

0 0 

66.0 


X 

PTC 

151 

85.3 

-64 0 

0 0 

0 0 

80.0 


X 

PTC 

152 

47.0 

0 0 

0.0 

0.0 

47 0 


X 

PTC 

153 

49.5 

28 6 

0 0 

0.0 

57 0 



PTC 

15A 

42.5 

46 5 

0,0 

0.0 

48,0 


X 

PTC 

155 

117.3 

-130.0 

0 0 

0.0 

101.0 


X 

PTC 

156 

132.4 

-145 3 

0 0 

0.0 

118 0 


X 

PTC 

157 

91 8 

-77.3 

0 0 

0 0 

86.0 


X 

PTC 

158 

70.0 

0 0 

0 0 

0 0 ' 

70 0 


X 

PTC 

159 

48.4 

16 2 

0.0 

0 0 

49 0 


X 

PTC 

163 

47 0 

0 0 

0 0 

0.0 

47 0 


X 

Cold Soak 

150 

80 0 

-3.9 

0,27 

0 0 

56 0 


X 










TABLE 3.3 (Continued) 


CURVE FIT COEFFICIENTS REPRESENTATIVE 
or THE CSS THERMAL ENVIRONMENT 


Mission 

Mode 

Node 

& 

o 

C"F) 


2 

(°F/Hr ) 

3 

<“F/Hr ) 

T 

inf 

Time 

Arjtument 

t+3 

(Hrs) 

■inH 


Cold Soak 

151 

115 2 

-IS 6 


-0 036 

31 0 


X 

Cold Soak 

152 

109.3 

-22 3 


-0 041 

5*0 


X 

Cold Soak 

153 

88.7 

-10 2 


-0 0074 

7 0 


X 

Cold Soak 

154 

107.5 

-18 3 


-0 03 

6 0 


X 

Cold Soak 

155 

60 7 

-0.25 

0.046 

-0.0021 

58.0 


X 

Cold Soak 

156 

51 3 

7 65 

-0 42 

0 0067 

96 0 


X 

Cold Soak 

157 

99.6 

-12.6 

1 07 

-0 029 

50 0 


X 

Cold Soak 

158 

50.0 

5.8 

-0.26 

0.0039 

96.0 


X 

Cold Soak 

159 

111.6 

-20 1 

1 42 

-0 035 

5.0 


X 

Cold Soak 

163 

32 1 

12 2 

-0 70 

-0.013 

109.0 


X 

Hoc Soak 

150 

103.3 

-10 2 

0.812 

-0.002 

50 0 


X 

Hot Soak 

151 

47.7 

8.73 

-0 443 

0 00071 

103 0 


X , 

Hot Soak 

152 

32 1 

12.2 

-0.7 

0.013 

109.0 


X 

Hot Soak 

153 

105,2 

75.8 

-882 4 

396 6 

105.0 


X 

Hot Soak 

154 

97 2 

62 1 

-870 3 

395 0 

99.0 


X 

Hot Soak 

155 

120.7 

40,4 

-975 4 

450 0 

122.0 


X 

Hot Soak 

156 

24 6 

17.9 

-0.892 

0.013 

136,0 


X 

Hot Soak 

157 

17.0 

17 6 

-1 2 

0 028 

105.0 


X 
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TABLE 3*3 (Continued) 


CURV^FIT COEFFICIENTS REPRCSENTATIVE 


Mleelon 
Mode , 

Node 

<’F) 

(°F/Ur) 

«2 

CF/ilr^} 

®3 

{'F/nr^) 

T 

inf 

°F 

Time ArRument 

t+3 

(Hts) 



Hot Soak 

158 

86 1 

-6.2 

0 2 

-0 0021 

25 0 


K 

Hot Soak 

159 

102.2 

63.1 

-685 9 

401 9 

10 0 


X 

Hot Soak 

163 

109 3 

-22.3 

1 7 

-0 041 

5.0 


X 

'lUnar Orbit 

150 

53 5 

5 5 

0 0 

0 0 

114 0 

X 


Lunar Orbit 

151 

74.6 

-1.52 

0 0 

0 0 

54 0 

\ 


Lunar Orbit 

152 

45 0 

0.0 

0 0 

0 0 

45 0 

\ 


Lunar Orbit 

153 

79.0 

0 0 

0.0 

0.0 

79 0 

* 


Lunar Orbit 

154 

53 0 

0.0 

0 0 

0.0 

58 0 

X 


Lunar Orbit 

155 

51 1 

2 63 

0 0 

0.0 

84 0 

X 


Lunar Orbit 

156 

53 9 

3 89 

0 0 

0 0 

107 0 

X 


Lunar Orbit 

157 

69,0 

0 0 

0 0 

0 0 

69 0 

X 


Lunar Orbit 

158 

75.72 

-2 24 

0 0 

0 0 

42.0 

X 


Lunar Orbit 

159 

57 0 

0 0 

0 0 

0 0 

57 0 

X 


Lunar Orbit 

163 

45 0 

0 0 

0 0 

0 0 

45 0 

■ 



*5 3 

r - K + a^X uhere X * time argument 














given by the following equations . 


■ “I - *»»cs - W <18a) 

<*r ■ (,\J - t/) (18b) 

g = AL3 (T^ - T^) + .9 (18c) 

Where: 

ALl = 01118x10“^ BTU/HR°r'^ 

AL2 = 1 25x10"^ BTU/HR°R^ 

AL3 = .37313x10"^ BTU/HR®R 

The value of AL2 will vary with actual tank heat leak characteristics. 

The amount and temperature of the hydrogen in storage is computed in the 
tankage model. Figure 3 13B presents typical hydrogen tank heat leaks 
for various flowrates and quantities. 

3.3 3.3 Vacuum Annulus Failure 

The nominal pressure in the tank annulus under the design vacuum 

-9 -7 

condition is approximately 2 x 10 psia (10 mm Hg) Depending upon the 
type of failure, the annulus pressure in the failed condition will be some 
where between the nominal pressure and the annulus burst disc pressure. 
Burst disc pressures for the oxygen and hydrogen tanks are approximately 
75 and 80 psia, respectively. Within these pressure bounds, gas in the 
annulus may exist in either rarefied or continuum regimes. The criteria 

for determining the gas regime is based primarily on the ratio of the 
molecular mean free path to the amount of space available for molecular 

motion This ratio is defined as the Knudsen number (K ) . For values of 

n 

0 01 the gas is normally treated as continuum flow, while for values 
of > 1.0 the gas is assumed to be in the rarefied regime. For the CSS 

_3 

tanks, an annulus pressure on the order of 10 psia represents the transi 
tion between the two regimes . 
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FIGURE 3 13A 


TYPICAL OXYGEN TANK HEAT LEAK 



FIGURE 3 13B 

TYPICAL HYDROGEN TANK HEAT LEAK 



TANK OUTFLOW ~ (LB/HR) 
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The type of analysis used to compute the amount of heat transferred 

* 

across the annulus depends upon the flow regime. For the rarefied regime 
the equation for determining the heat flux per unit area is given by 
Reference 5 as: 


q - 



(19) 


For continuum flow regimes, the heat flux per unit area is given by Reference 
6 as: 




( 20 ) 


An analysis was performed for the Apollo cryogenic storage tanks 
(Reference 7) assuming residual gas in the annulus consisted of air or of the 
appropriate cryogen. The analysis showed that for an annulus pressure of 
less than approximately 10 ^ psia, heat conduction by the residual gas was 
essentially negligible. However, heat conduction increased rapidly with 
Increasing annulus pressure. The total heat transfer to the tanks as a 
function of the annulus pressure is shown in graphical form in Figures 3.14 
and 3.15. These curves have been incorporated into the thermal model such 
that when'the thermal subroutine encounters a failed tank flag, the data 
presented in the figures is used to determine the heat flux to the tank based 
upon the vacuum jacket pressure. This pressure is an input to the program. 

3, 3. 1.4 Components and Line Environmental Temperatures 

The enviornmental temperatures used to determine the state of the fluid 
in the components and lines were obtained in a manner identical to the tank 
environmental temperatures. The thermal mathematical model of the Apollo 
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TOTAL TANK HEAT LEAK (BTU/HR) 



- 10.0 - 8*0 - 6.0 - 4.0 - 2.0 0 2.0 

Logj^Q P, Annulus Pressure (Psia) 
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TOTAL TANK HEAT LEAK (BTU/HR) 


FIGURE 3.15 - VARIATION OF CSS TANK HEAT LEAK FOR LOSS OF 
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Service Module was used to obtain data at several locations as gxven in 
Table 3.4. These data were incorporated into the Integrated Systems 

4 

Program as a series of curve fit correlations. The form of the resulting 
equation is the same as Equation (18). The coefficients for the curve fit 
equations are presented in Table 3.3. 

The TMM did not model the Command Module (CM) cabin, the CM Reaction 
Control System Tank Annulus, or the ECS umbilical. The temperature at these 
locations was assumed to be 70®F. 

3.3.2 Thermophysical Properties for Oxygen and Hydrogen 

Computer subroutines for computing the thermophysical properties for 
oxygen and hydrogen were developed and included into the Integrated Systems 
Program. The majority of the data used in these subroutines was recommended 
by the National Bureau of Standards (NBS). That data not available from 
the NBS was taken from Reference 8 (Stewart) . In general these subroutines 
determine pressure ~(p) , temperature (T) , density (p), enthalpy (h) , the 
specific heat at constant pressure (c^) , the specific heat at constant 
volume (c^) j viscosity (y), thermal conductivity (k) , isotherm and isochor 
derivatives and the functions (j> and 0, providing sufficient information 
concerning the fluid state is known. The data is included in the subroutines 
both in equation and table form. A routine is also included in the subroutines 
for interpolating between values from the tables. 

3. 3. 2.1 Oxygen Properties 

Data contained in the subroutines for computing the thermophysical 
properties of oxygen were taken from References 8, 9, and 10. Reference 8 
presents an equation of state for oxygen determined by R. B. Stewart. The 
Isotherm and isochor derivatives were derived from this equation. While 
Reference 9 (Weber) also presents data for the isotherm and isochor 
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TABLE 3,4 


TEMPERATURE INDICATOR LOCATION 


SINK INDEX 

LOCATION 

150 

TUNNEL 

151 

Bay 4 

152 

Bay 4 

153 

FoTOard Bulkhead 

^ 154 

Bay 4 

155 

Bay 4 

156 

Bay 4 

157 

Bay 4 

158 

Forward Bulkhead 

159 

Forward Bulkhead 

160 

Umbilical 

161 

CM RCS Tank Annulus 

162 

CM Cabin 

163 

Bay 1 
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derivatives, the results derived from Reference 8 are valid for a greater 
range. Both sets of data were included in the thermophysical subroutines 
to provide greater flexibility since the derivatives derived from Reference 
8 may be determined with knoTO values of temperature and density, while 
the data from Reference 9 uses either pressure and temperature or pressure 
and density. 

The equation of state provided by Reference 8 may be used to compute 
pressure from known values of density and temperature. This equation is 
valid for temperatures between 117 and 540 “R and for pressures to 5000 psia. 
The equation of state is as follows: 

p = pRT + (n^T + n^+ n^/T^ + n^/T^ + n^/T^)p^ 

+ (ngT^ + n^T + ng + n^/T + n^g/T^)p^ 

+ p (n^g/T + n^g/T + " 20 ^^ ^ ^^25p ^ 

+ p^(n2^/T^ + ^^25p^^ 


^^ 28 ’^^^ , ^28 


‘28. . . "28 
)exp[h2g(P -P 


( 21 ) 


where: 


= 13.333 gmol/liter, critical density 
= 154.77°K, critical temperature 


In this expression, p is in atmospheres, T is in ®K, and p is in gmol/liter. 
The constants for this equation are given in Table 3.5. The subroutine in 
which Equation (21) is programmed converts the units of density from Ibm/cu-ft 
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CONSTANTS FOR 



E 3.5 


'S EQUATION OF STATE 

1.73655508 x 10^ 

3.01752841 x 10^ 

n^^=-3. 49528517 x 10^ 

n^g= 8.86724004 x 10~ 

n^g=-2. 6781766 7 x 10^ 

H 2 q= 1.05670904 x 10^ 

5.63771075 x lO" 

n22=-l- 12012813 

ii 23 = 1.46829491 x 10^ 

9.98868924 x 10“' 

n23=”0. 00560 

n„ =“0.157 
26 

n2y=-0. 350 
ii2g= 0.90 



to gmoi/liter and temperature from “R to °K prior to use. The output 
pressure is xn psia. 


The xsotherm derxvatxve, 


i£| 

3p' 


, xs the partxal derivative of Equatxon 


(21) wxth respect to the density at constant temperature and is given by! 


ap 


24 

i = R T + I n X 
1=1 ^ ^ 


( 22 ) 


Values for the variables in this equation are presented in Table 3.6. 

The coefficients n^ are the same as those for Equation (21) . The 
subroutine which computes the isotherm derivative converts density from 
Ibm/cu-ft to gmol/liter and temperature from “R to ®K prior to its use in 

3 

Equation (22). The output derivative is in units of psia-ft /Ibm. 

The isochor derivative |^| ^ , is the partial derivative of Equation (21) 
with respect to temperature at constant density and is given by: 


i£| 


ST'p 


pR + 


24 

I 

1=1 


n X 

1 X 


(23) 


where the coefficients n.^ are the same as those listed for Equation (21). 

The variables in this equation are shown in Table 3.7. The conversion from 
one set of units to the other is performed the same as for Equations (21) 
and (22) . The output derivative is in units of psia/^R. 

In addition to providing data for the isotherm derivative and the 
isochor derivative. Reference 9 presents data for density, enthalpy, 
specific heats (at constant pressure and volume) and the functions 0 and 
4>. All are in tabular form. Further flexibility was provided by permitting 
these properties to be determined from known values of pressure and temperature 
or from pressure and density. These tables are valid for temperatures from 
100 to 540°R at 5 degree increments and for pressures ranging from 367 to 
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TABLE 3.6 

VARIABLES FOR COMPUTING ISOTHERM DERIVATIVES 



= 2pT 

^13 

= Sp'' 

^2 

= 2p 

^14 

= 5p^/I 

^3 

= 2p/T^ 

^15 

= ty 

^4 

= 2p/T^ 

^16 

- 

^5 

= 2p/T^ 

^17 

- 


= 3pV 

^18 

« fy/T^ 

^7 

= 3p^T 

^19 

- fy/T^ 

^8 

= 3p2 

^20 

= f^/T^ 

^9 

= 3p^/T 

^21 

= fg/T^ 

\o 

= 3P^/T^ 

^22 

= fg/T^ 

\l 

3 

= 4p T 

^23 

- ^8/^' 

^12 

^24 

- 

= fgf2^3 ■'■ 

^”28+l' , . /"28+l' 

p + P fjfy 

exp 

(n25 P ) 

*6 

94= 2 
« 3f^p ■ 

p"23 _p "28 

^9 

h 

= 5f^p 


= exp[n2gf2 + n^/T-T^)^] 


2n2?(T-T^)f2 


= 2f3^pn25 


fg = 7f^p^ + f3p^ 

^9 " ^^28+1^ P 
^10 " 

^11 ^^ 2 ^ 26 ^ 3^10 
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TABLE 3.7 


VARIABLES FOR COMPUTING ISOCHOR DERIVATIVE 




2 






■= P 

" ^13 


0 


^2 

= 0 

^14 

= 

5/„2 
-p /T 


^3 

= -2p^/T^ 

^15 

= 

-2p'fl/T^ . 



= -4p^/T^ 

^16 

= 

-Sp^f^/T^^ 


X 

5 

= -6P^/T^ 

^17 

= 

-4pf^/T^ 


^6 

3 

= 2Tp 

^18 

= 

-2p^f^/T^ 


X 

3 

- P 

^19 


-3p^ 


^8 

= 0 

^20 

= 

-4p^f^/T^ 



= -p^/T^ 

^21 

- 

-2p^f^/T^ 


^10 

- -2p^/T^ 

^22 

= 

-3p'f,/T< 


^11 

4 

= P 

^23 

= 

-4p^fj^/T^ 



^24 


-^’^28+1^ , 


^12 

= 0 


0 Va 


where: = exp (n 2 gp 

. 





. ^28 

^28 





tj -p -p^ 
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1250 psia at increments of 73.5 psi. Table 3.8 lists the thermophysical 
parameters which are 'available, 

Bata from Reference 10 is also in the form of tables and provides 
viscosity and thermal conductivity as functions of pressure and temperature. 

This data is valid for temperatures ranging from 180 to 540°R at increments 
of 18 degrees and for pressures from 367 to 1250 psia at increments of 73.5 psi. 
3. 3. 2. 2 Hydrogen Properties 

The thermophysical properties used in the hydrogen subroutine are 

based upon computer programs developed by the NBS (Reference 11) . These 

properties are valid for pressures from 1 to 5000 psia and for temperatures 

ranging from 25 to 5000“R. In the NBS programs, all thermophysical properties 

are determined from known values of pressure and temperature and include 

- enthalpy, density, specific heats for both constant pressure and constant 

volume, viscosity and thermal conductivity. However, in the subroutine 

developed for the Integrated System Program, provision was made to compute 

the pressure and temperature from other known properties as described below. 

3 

The pressure, density, and temperature are in units of psia, Ibm/ft , and 
“R, respectively. 

The HBS program expresses hydrogen density as a polynominal equation 
in terms of pressure and temperature, p = p(p,T). The equation has been 
programmed such that pressure may be determined from this equation as a 
function of density and temperature using a Newton-Raphson iteration method. 

An a priori value of pressure is required to initiate the iteration and 
should be estimated as closely as possible to the actual value. Temperature 
is computed in a manner similar to that of pressure from this same expression. 
The available parameters are listed in Table 3.8. 
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TABLE 3.8 


THERMOPHYSICAL PROPERTIES SUBROUTINES 


- PROPERTY 

UNITS 

INPUT 

VARIABLES 

REFERENCE 

OXYGEN 

Density 

Ibm/ft^ 

P&T 

9 

Enthalpy 

btu/lbm 

P&T 

9 

Enthalpy 

btu/lbm 

P&p 

9 

Function 0 

btu/lbm 

P&T 

9 

Function (j> 

psia-ft /btu 

P&T 

9 

Isochor 

psia/“R 

T&p 

8 

Isochor 

psia/°R 

P&T 

9 

Isochor 

psia/°R 

P&p 

9 

Isotherm 

3 

psia-ft /btu 

T&p 

8 

Isotherm 

psia-ft^/btu 

P&T 

9 

Iso therm 

3 

psla~ft /btu 

P&p 

9 

Pressure 

Psia 

T&p 

8 

Specific Heat, 

s 

btu/lbm- 

P&T 

9 

Specif xc Heat, 

S 

btu/lbm-®R 

P&p 

9 

Specific Heat, 
C 

V 

btu/lbm- °R 

P&T 

9 
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TABLE 3.8(coNT'D) 

THERMOPHYSICAL PROPERTIES SUBROUTINES 


PROPERTY 

UNITS 

INPUT 

VARIABLES 

REFERENCE 

OXYGEN 

Specific Heat, 

btu/lbm-°R 

P&p 

9 

C 




V 




Therml Conduc- 

btu (ft/hr)/ 

P&T 

10 

tlvity 

ft^/“R 



Temperature 

“R 

P&p 

9 

Viscosity 

Ibf-hr/ft^ 

P&T 

10 

HYDROGEN ] 

Density 

Ibm/ft^ 

P&T 

11 

Enthalpy 

btu/lbm 

P&T 

11 

Function 0 

btu/lbm 

P&T 

11 

Function <J) 

3 

psia-ft /btu 

P&T 

11 

Isochor 

psia/®R 

P&T 

11 

Isotherm 

3 

psia-ft /btu 

P&T 

11 

Pressure 

Psla 

T&p 

11 

Specific Heat, 

btu/lbm-“R 

P&T 

11 

C 




P 




Specific Heat, 

btu/lbm-°R 

P&T 

11 

C 




V 




Thermal Conduc- 

btu(ft/hr)/ 

P&T 

11 

tivlty 

ft^/“R 



Temperature 

“R 

P&p 

11 

Viscosity 

Ibf-hr/ft^ 

P&T 

11 
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4.0 CONCLUSIONS 


The Apollo Cryogenics Integrated Systems Program was developed to aid 
in the detailed study of preflxght analysis and postflight evaluation of 
the Apollo Cryogenics Storage System. Capable of accurately simulating 
system performance, the program is very flexible, particularly in the area 
of malfunction analysis. All significant effects on system operation have 
been included. The program will greatly enhance the ability to analyze 
the performance of the Apollo Cryogenic Storage System. 

The program has been used for the Apollo 14 preflight analysis of the 
CSS performance. Simulations of the nominal mission, system operation 
with an anomalous oxygen check valve, and system performance with several 
-cases of degraded vacuum in the annulus of the oxygen tank were conducted. 
These analyses, however, were carried out before the effects of oxygen 
thermodynamic stratification were incorporated into the program and are 
therefore not presented in this document. 

It is planned that the program will be used in the preflight analysis 
of the Apollo 15 Mission. The effects of thermodynamic stratification of 
oxygen will be considered for this study. In addition, prediction of 
heater temperatures will be included. Reports will be issued which docu- 
ment both the nominal mission performance and the results of several mal- 
function cases. 
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APPENDIX I 
NOMENCLATURE 


Area , in , 

Accommodation coefficient, equal to 1.0. 

V h , 1 

Capacitance analog, C - , Ibm - in /Ibf 

Specific heat of fluid at constant pressure, btu/ 
(lbra-°R). 

Specific heat of fluid at constant volume, btu/ 
(lbm-®R). 

Wall thickness, in. 

Line diameter, in. 

2 2 4 

Electrical potential analog, E = p , Ibf /in 
Energy of system 

Friction factor, dimensionless. 

2 4 

Conductance analog Ibf -hr/(lbm-in ) 

2 

Gravity constant Ibm-in/ (Ibf-hr ) 

Enthalpy of stored fluid, btu/lbia. 

Enthalpy at initial state, btu/lbm 

Enthalpy of fluid in vapor cooled shroud, btu/lbm. 

Current analog, I = m, Ibm/hr. 

Residual current analog, I_ = Ibm/hr. 

^ RT^ 

Thermal conductivity, btu/ (hr-ft-“R) . 

/g n>.l/2, ..v n+1 ^ 

Constant, K =(— \~T/ 2(n-l), lbm-°R /(Ibf-hr), 


Knudsen number, X/1, dimensionless. 
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1 

L 

M 

m 

m 

n 

N 


V 

q 

q 

Qr 

r 

R 

s 

t 

T 


vcs 


Vi 

u 


V 

V 


Characteristic length, ft. 

Length, in. 

Gas molecular weight, Ibm/ (lb-mole) 

Mass, Ibm. 

Mass flowrate, Ibm/hr. 

Ratio of specific heats, c„/c , dimensionless. 

P V 

Poly tropic exponent for gas lines attached to tank, 
dimensionless . 

Reynolds number, dimensionless. 

Pressure, psia 

Total heat input rate, btu/hr. 

Heat leak flux, btu/ (hr-ft ) 

Radiation heat leak, btu/hr. 

Mean pressure vessel radius, in. 

Universal gas constant. 

Inner tank wall radius, ft. 

Outer tank wall radius, ft. 

Fluid velocity, in/hr. 

Time , hours 
Temperature, °R, 

Temperature of tank outer wall, "R 
Temperature of tank inner wall, °R. 

Average temperature of fluid in vapor vooled shroud, "R. 
Internal energy of fluid, btu/lbm. 

2 

Overall heat transfer coefficient, btu/ (hr- In -®R) . 

3 

Specific volume, ft /Ibm. 

3 

Tank Volume, ft . 
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Y 
a 
3 
A 
£ 

6 

X 

P 

P 

a 

V 

Subscripts 

b 

c 

h 

1 

3 

k 

1 

LM 

r 

R 

s 

T 


Volume of lines attached to tank, ft . 

Young's modulus for tank material, psi. 

Tank material thermal expansion coefficient, 1/®R. 

Tank thickness. In. 

Difference, dimensionless 

Emissivity, dimensionless 

Thermodynamic function, , btu/lbm. 

dp I p 

Mean free path, ft 
Viscosity, Ibm/ (in-hr) . 

3 

Fluid density, Ibm/ft 

-11 

Stephan-Boltzman constant, = 1.1x10 , btu/ (hr-in 

°R^). 

Poisson's ratio for tank material, dimensionless. 

i iE 

Thermodynamic function, p 9u 
Line bend angle, degrees. 

Bulk fluid 
Cracking 
Heater tube 
Inlet 
Outer 

Arbitary node or line 
Arbitary node or line 
Logarithmic Mean 
Reseat 
Residual 

/ 

Sink 
Thermal 
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APPENDIX II 


DERIVATION OF dp/dt FOR VARIABLE VOLUME TANKS 
DUE TO PRESSURE AND TEMPERATURE CHANGES 


Introduction . 

An expression for the rate of change of pressure for the cryo- 
genic storage tanks has been developed for a constant volume 
(Reference A-1) and for a variable volume caused by tank stretch due to 
pressure changes (Reference A-2). The latter expression is; 


where ; 


iE 

dt 



1 + 30|pr 
2dY 


q 

(l-u) 


0 


= - 


P 

P 


(AII-1) 


(AII-2) 

(AII-3) 


If Young's Modulus, Y, is taken as infinity (no tank stretch). Equation 
(AII-1) reduces to the expression for dp/dt as -presented in Reference 
A-1. 


The expression for dp/dt considering variable volume resulting 
from tank stretch due to both pressure and temperature changes is 
developed in this Appendix. For completeness, a portion of the 
derivation from Reference A-2 is also presented. 


Derivation ♦ 

The First Law of Thermodynamics may be written: 


dEn 

dt 


q + m h' 



If It IS assumed that the system is at uniform pressure, p, 
and that kinetic and potential energy may be neglected, then 


En — mu = pVu 
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and 

Now 

Then 


Since ha: h' 

Thus , ^ = 
dt 

Assume that 
then 


d(f>Vu) 

dt 


q + m h' 



udfoV) + - q + ; h’ - 


u = h - pv. 

udm . oVd(h-pv) . ' , i pmdv pvdm 

“d? - q + - h’ - ^ ^ 

• * ^ ♦ 

um + pV - pV d (pv) = q + m h'- pm dv - p v m 
dt dt dt 

* * dv 

m (u + pv) + pv ^ - pVp ^ - pVv ^ + ra h' - pm — 

dt dt dt 


m h + pV ^ - m ji^'av - V ^ = q 4- m h' - 
dt dt dt 


V^=pV^ -q+m(h-h') 
dt dt 

dp = p ^ - q_ + m 
dt dt V V (h - h'). 

, the last term is negligible. 



>Qh ~ d 
dt V 


(AII-4) 


h = h(p,p). 


dh 


Ih 

8p 


dp + ^h_ 
P 8p 


P 


dp 
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= 8h 


+ 8h 

dp 

3p 

p dt 

9p 

p dt . 


How, 


and 



d f— I 

dp ^ JV ^ 1 ^ 

dt dt V dt 


ni dy _ £. dm _g^ dV 

^2 dt m dt V dt . 


So, 


dt 


% 

3p 


d£ 

+ 3h 

p dm ' 

- p dV‘ 

P dt 

ap 

1 p _ m dt 

V dt . 


Substitutxng for0 from Equation (AII-2) 



dh 

_ 3h| 

1 dp - 

^ ra + 


dt 

9p 1 p dt 

m V 

Since 

h^= 

u + 

pv = 

= u 

P , 

then 



_9u 

p + 

3(”) 
1 


8p 

c= 

p 

ap 

3p Ip 


3h : 


8u 

. 1 . 

1 


3p i 

p 

3p 

•r 

P 

p 



Substituting for G from 


dt 



1 1 dp 0 ' 0 dV 

p J dt m ^ V dt . 

Equation (AII-3) 

i] iE _ ® ® iZ 

p J dt ra ^ V dt . 


(AII-5) 
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It xs at this point that the assumption is made that tank stretch 
is a function of both pressure, p, and temperature, T. 


^ ^ iv U + iv 

dt 3p dt 9T 


p dt 


(AII-6) 


It Is necessary to determine dT/dt 

dT _ ^1 dp 3T I dp 

dt 3p Ip dt 9p ip dt ; 


(All- 7) 


Now, m - PV 


dm „ dp ,, dV 
dt dt ^ dt 


So, dp _ m P dV 
dt " V " V dt 


(AII-8) 


Substituting Equations (AII-8) into Equation (All-?), and substituting 
Equation (AII-7) into Equation (AII-6) results in* 

^ = 3V I ^ r^l ^ 

dt 9p It dt BTlpLsplp dt 9plp V ~ V 

For convenience, the following substitutions have been made: 

9V 


9p 


B = 


3T 


C = 


9T 

9p 


Then, 


D = 


3p 


~ = A^ + B Tc^ + D ^)1 

dt ^dt ^ r dt ^ V dt-'J 


~ = (A + BC) ^ + BD - - ^ — 
dt ^ dt ^ V V dt 
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dV 

dt 


= (A + BC) 


(1 + 


BDP 

V 


iE + ^ 

dt V 


dt 


(A + BC) 


1 + 


dt 

BDP 

V 


BBm 

V 


(AII-9) 


Substituting Equation (AII-9) into Equation (AII-5). 

(A + BC) 


dt 


1 

p(j) p 


^ _ 

dt 


0 


m 


m 


G 

V 


dt 




1 + 


BPp 

V 


] 


it _ r 1_ 4. i + Q(A + BO ] ^ ^ q. QBDct 

dt [ pij) p V+BDp J dt m V(V+BDp) . 


But from Equation (AII-4) 

dh 1 dp , q 
d? "p ^■^pV ■ 


Therefore, 

• • 

, l/ e (A+BC) ~| d£^ 0m 0BDm 

(/dt pV “[ p<i) / V+BDp J dt m V(V+BDp) 

and, 

• « 

r 1 e(A+BC) 1 dp , 0m 0EDm ^ ^ 

L p<j> V+BDp J dt m V (V+BDp) pV 

resulting in: 


0tn 

dp ^ pV 
dt 


1 

"p? 


■ q _ 0BPm 
PV V(V+BDP) 


+ e . (A+B .C) 
V+BDp 
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or 


Oii({ ) , _ 0(j)BPnip 

i£ = V V V(V+BDo) 

1 + 600 (A+BC) 

V+BDp 


Substituting for A, B, C and D 


eM 

V 


iR = 

dt 


1 + 


+ _ 
V 


04*mP 


3V 


3T 


V(V+ 


TT lp~3^ 

av 1 St I 


p St 


Spli 




V + p- 


T ^ St 


p 9p 


) 


sv 


ST 


ST Ip Sp Ip 


(AII-10) 


Reference A-2 presents the derivation for change in volume with respect 
to pressure at constant temperature as. 



3rV 

2dY 


(1-u) 


(AII-11) 


The change in volume with respect to temperature may be obtained as 
follows : 


iy. 

St 

jv 

Sr 


1 sv 

3r 





I 

a>l 

1 

P ST 1 

P 




a 

m3) 

/ 2 

= ly 


3r 

= 4 irr 

r 

* 

s in temperature. 

then 




ST 


r AT 


= 3dV 


(AII-12) 


74 



Substituting these expressions into Equation (All- 10) results ml 


V V 


Oijiiii 3 ctp -r— 

^ 


= 

dt 


V (l+p3ot3p 


Ip) 


1 + 


04.P(|S^ (1 


u) +3o|^ 
op 


1 + 3ap ^ 


(AII-13) 


It IS necessary to account for one additional effect which results in 
a variation in the change in volume with respect to pressure. This 
effect IS the compression of the fluid in the fill and vent lines of 
the tank. These lines extend from the pressure vessel to the warm 
environment. The fluid at the cold ends is relatively incompressible 
as compared to the near ideal gas conditions at the other ends. When 
pressure increases , the fluid in the lines is compressed by the 
expanding stored fluid which move up the lines. To account for this 
effect, it was assumed that the lines will be filled with an ideal gas 
at ambient temperature and that no mixing between the tanked cryogen 
and the perfect gas would take place. Then 


Thus, 


or 


The mass flow. 


pV^ = constant. 

Li 


iE 




n-l dV. 


L 0 


dt 


dVi 

dt 


-V. 


k iE 

np dt 


m, is the total flow from the tank: 


(All- 14) 


m = m + m., (AII-15) 

o J. 

where: m^ = flow to EPS and ECS subsystems, Ibm/hr 

• dV 

= p L = flow from fill and vent lines, Ibm/hr 

dt 
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Substituting these expressions into Equation (AII-10) results in: 


6m({> 

V 

ik = 

dt 

1 + 


<j>q 

V “ 


G([im 3ap 
oP p 


M. 

V (l+p3a9p 


Ip) 


Q<{>P(||y (1 - u) +3ag[p ) 


1 + 3ap If 

9p p 


(AII-13) 


It IS necessary to account for one additional effect which results in 
a variation in the change in volume with respect to pressure. This 
effect IS the compression of the fluid' in the fill and vent lines of 
the tank. These lines extend from the pressure vessel to the warm 
environment. The fluid at the cold ends is relatively incompressible 
as compared to the near ideal gas conditions at the other ends. When 
pressure increases, the fluid in the lines is compressed by the 
expanding stored fluid which move up the lines. To account for this 
effect. It was assumed that the lines will be filled with an ideal gas 
at ambient temperature and that no mixing between the tanked cryogen 
and the perfect gas would take place. Then* 


Thus, 


pV = constant. 

Lr 


V. 


n dp 
L dt 


+ npV! 


n-1 dVj^ 
dt 


= 0 


or 


The mass flow. 



dt np dt 

« 

m, is the total flow from the tank; 


(All- 14) 


m = m + m„ (AII-15) 

o 2. 

where: = flow to EPS and ECS subsystems, Ibm/hr 

dV, 

m = p 1 = flow from fill and vent lines, Ibm/hr 

dt 
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Note* Flow out of the tank is negative in sign 
Substituting Equation (AII-14) into Equation (AII-15) yields: 

L dp 

m = m — p — 

o np dt 


l>Jhen this equation is substituted into equation (AII-13) and rearranged, 
the resulting expression is: 


4- ‘i’ 
V 


0ip m 3ap— 
— _o 9p_ 


(1 + p 3a 


3T 

9p 


dt 


T , 0(j) pV, 0<j) 3ap" 

np 


9T 

9p 


where : 


np (1 + 3ap 

0(j.p(^Y Ip^ 


9p 


+ F 

) 


1 + 3ap 


Sp 


(AII-16) 


Thus, Equation (AII-16) describes the rate of change of pressure for 
a variable tank volume considering both pressure and temperature 
changes. Note that if the coefficient of thermal expansion, a, is 
equal to zero and if the plumbing fluid compression volume, Vj^, is 
zero, Equation (AII-1) is obtained. Table A-1 presents the 
constants used for both the oxygen and hydrogen tanks . The partial 
derivatives are obtained by differentiating the equation of state. 
Since the Integrated Systems Program uses an equation of state that 
is a function of temperature and density, it is desirable to obtain 
8T/9p in the following form. 

^ 9p I _ 1 

8p Ip 3p T 3T Ip 

This relation was taken from Reference A-4, page 32, Equation (3-11). 
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Thus, this partial derivative may also be obtained from the equation 
of state. 
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TABLE A-1 

CRYOGENIC TANK PROPERTIES 


OXYGEN 

Material. ... Inconel 718 

Volume, V, ft 4.75 

Inside Radius, r, in 12.513 

Young's Modulus, Y, psi 30 x 10 

Poisson's Ratio, u, dimensionless 0.29 

Wall Thickness, d, in .....0.059 

Coefficient of Thermal Expansion, a, /®R 5.93 x 10 ^ 


HYDROGEN 

Material 

Volume, V, ft^ 

Inside Radius , r, in. 

Young's Modulus, Y, psi 

Poisson's Ratio, u, dimensionless 

Wall Thickness, d, in 

Coefficient of Thermal Expansion, a, /°R 


Titanium-5AL-2.5SN ELI 
6.80 
14.103 
17 X 10^ 

0.30 

0,044 

1.23 X 10“^ 


Note: Values for a were derived from Reference A-3, Table 7.3.4, 

page 7-59. 
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TABLE A-1 

CRYOGENIC TANK PROPERTIES 


OXYGEN 

Material Inconel 718 

Volume, V, ft^ A. 75 

Inside Radius, r, in 12.513 

Young's Modulus, Y, psi 30 x 10 

Poisson's Ratio, n, dimensionless..... 0.29 

Wall Thickness, d, in 0.059, 

Coefficient of Thermal Expansion, a, /“R 5.93 x 10 ^ 


HYDROGEN 

Material 

Volume, V, ft^. . . 

Inside Radius, r, in 

Young's Modulus, Y, psi 

Poisson’s Ratio, u, dimensionless 

Wall Thickness, d, in 

Coefficient of Thermal Expansion, a, /“R 


Titanium-5AL-2.5SN ELI 
6.80 
14.103 
17 X 10^ 

0.30 

0.044 

1.23 X lO"^ 


Note; ■ Values for a were derived from Reference A-3, Table 7.3.4, 
page 7“59. 
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APPENDIX III 


DERIVATION OF EQUATIONS USED IN THE 
COMPONENT AND LINE MODEL 


MATRIX MODEL 

The matrix model determines the pressures for each of the oxygen and 
hydrogen system nodes shown, in Figures 3.5 and 3.6 respectively. 

The conservation of mass principle may be written; 

^ m = . (AIII-1) 

For a perfect gas, 

p = pRT. CAIII-2) 


Combining Equations (AIII-1) and (AIII-2) and performing the 
differentiation results in: 


m 


(T* - pg) . 


RT 


2 ' dt ^dt 


To set up an electrical analogy, let 

E E 

for an equivalent driving function, and 

I = m 

for the current. 


(AIII-3) 


(AIII-4) 


(All 1-5) 


In terms of the transforms variables, Equation (AIII-3) becomes: 


or 



V 

2RTv^ 


dt ■ dt 



R 


where 


C 


V 

2RT>^ 


(AIII-6) 


(AIII-7) 


(AIII-8) ' 
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(AIII-9) 


IS an equivalent capacitance, and 

T = ^ 

R rt2 

is an equivalent residual current . 


NODAL MODEL 

-The nodal model determines the pressures at selected locations 
along each leg of the oxygen and hydrogen systems. Figures 3.7 and 
3.8 respectively. 


The instantaneous pressure drop for incompressible flow through a 
line with inlet i and outlet j is given by the Darcy-Weisbach equation. 


p. - p 

1 3 


ps^ 

^ ij 2S, 


For incompressible flow in a circular tube, 

trD^ps 

TO = pAs = — . 


(AIII-10) 


(AIIl-11) 


Combining Equations (AIII-10) and (AIII-11) and simplifying results in 

8 m^ 

'3 ^13 


p - p. = (tt-' 


(AIII-12) 


Let the mean density between 1 and 3 be based on the average pressure 
and temperature. Then, 

1 P^+P 

^13 E. 31 .j-y (AIII— 13) 

3. 3 • 


Substituting Equation (AIII-13) into Equation (AIII-12) and simplifying 
yields: 



8R(T +T.) 

-i 3 

g^r^D^ 


'■ D ' 




ij 


For laminar, incompressible flow in a circular tube. 


(AIII-14) 


f = 


16 


4:rpD 


(AIII-15) 


RE m 

Combining Equations (AIII-14) and (AIII-15) and simplifying results in: 
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X28Xf|iR ^ 

g TTD^t 1 J 


m 


(AIII-16) 


In terms o£ the electrical analog variables, Equation (AIII— 16) becomes: 


where 


E 


i 


E 

J 



(All 1-17) 


G = 


128LyR 

g^irD^ ^ 1 


(AIII-18) 


An Instantaneous heat balance through a line with inlet i and out- 
let j is: 


UA = m c (T -T ) . 

XjM p j ^ 


(AIII-19) 


where 


AT. 


LM 


(T^_T^) - (T^-T^) 


In 


’T, -T ' 
k 1 

T, -T 


(AIII-20) 


Combining Equations (AIII-19) and (AIII-20) and simplifying yields: 



J 


(AIII-21) 
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